Markovian and Post-Markovian dynamics of nonequilibrium thermal entanglement 
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The dynamics of two spins coupled to bosonic baths at diflterent temperatures is studied. The an- 
alytical solution for the reduced density matrix of the system in the Markovian and Post-Markovian 
case with exponential memory kernel is found. The dynamics and temperature dependence of spin- 
spin entanglement is analyzed. 



I. INTRODUCTION 

The influence of the environment plays an essential role in the description of the realistic quantum system 
Usually, environment destroy entanglement in the subsystem of interest. However in some cases it can create quantum 
correlations in the system 0, H, 01 • One of the ways to understand the role of the parameters of the system is to study 
exactly solvable models. Here, we study the dynamics of a model that was recently introduced by L. Quiroga 5]. It 
consists of two interacting spins in contact with two reservoirs at different temperatures. In such a non-equilibrium 
case most studies are restricted to the steady-state solutions [1, 0, 0, Q or to the zero temperature limit Q . 

This paper is organized as follows. In Sec. II we describe the model of a spin chain coupled to bosonic baths 
at different temperatures and derive a master equation in Born-Markov approximation. In Sec. HI we present the 
analytical solution for the system dynamics in the Markovian case, details of the solution can be found in Ref. [loj . 
In Sec. IV we present the analytical solution of the Post-Markovian master equation recently introduced by Shabani 
and Lidar Finally, in Sec. V we discuss the results and conclude. 



II. MODEL 



We consider a system of two interacting spins, with each spin coupled to a separate bosonic bath. The total 
Hamiltonian is given by 

H — Hs + Hbi + Hb2 + HsBl + HsB2, 

where 

Hs = |<Tf + |<7| + Kia+a^ + a^a+) 

is the Hamiltonian describing spin-to-spin interactions and &^,af are the Pauli matrices. Note, that the units are 
chosen such that kB = h = 1. The constants ei and £2 denote the energy of spins 1 and 2, respectively and K denotes 
the strength of the spin-spin interaction. The Hamiltonians of the bosonic "baths" for each spin j = 1,2 are given by 

n 

The interaction between the spin subsystem and the reservoir with creation operators ■ is described by 

n n fj. 
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The operators of the transitions in dynamical subsystem V,_^ are chosen to satisfy [Hs, Vj,ti] = '^i./j'^i./i; ^-nd the /^^^ 
act on the reservoir degrees of freedom. The total system (two spins with reservoirs) is described by the Liouville 
equation 

We assume that the evolution of the dynamical subsystem (coupled spins) does not influence the state of the envi- 
ronment (bosonic reservoirs) so that the density operator of the whole system a{t) can be written as: 

d(t) =p(i)Bi(0)B2(0), 

where each bosonic bath is described by a canonical density matrix Bj — e^'^^^^J /tr[e^^^^^^] and p{t) denotes the 
reduced density matrix of the spin subsystem. 

In Born-Markov approximation the equation for the evolution of the reduced density matrix [l^ is: 

, with dissipators 

and where the spectral density is given by 



JO 

To find a solution we go to the basis of the eigenvectors \Xi) with eigenvalues Xi of the Hamiltonian Hs, 

|Ai) = |0,0), Ai = -^i±^, 



IA3) =cos(0/2)|l,O) +sin(0/2)lO,l), A3 = k, 



IA4) = -sin(6i/2)|l, 0) + cos(6l/2)|0, 1), A4 = -k, 



where k = ^ K"^ + and tan6' = 2K/{Ae). In this representation the dissipators Ci{p) becomes 
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with transition frequencies 

uii = A2 — A3, 



■^2 — A2 + A3 
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and transition operators 

1/1,1 = cos(^/2)(|Ai)(A3| + |A4)(A2|), 
1/1,2 = sin(^/2)(|A3)(A2|-|Ai)(A4|), 
1/2,1 =sin(0/2)(|Ai)(A3|-|A4)(A2|), 

V2,2=COs(^/2)(|A3)(A2| + |Ai)(A4|). 

In this paper wc consider the bosonic bath as an infinite set of harmonic oscillators, so the spectral density has the 
form J^^\uj^,) = 7j(cj^)nj(a;^), where nj{uj^) = (e^^'^f - and J^^\-ojf_,) = e'^^'^f J^^\ojf_,). For simplicity we 
choose the coupling constant to be frequency independent ji{uj) = 7i and 72(w) = 72. In the basis |Ai) the equation 
for the diagonal elements of the reduced density matrix is given by 



It 



(Pii{t)\ /Pii(t) 

P22{t) \ \ P22{t) 

P33{t) I I P33{t) 

\ p4A{t) j \ pii{t) 



where i? is a 4 x 4 matrix with constant coefficients. The time-dependence for the non-diagonal elements has the 
following form 

Pi,,(t)=e*^-^p,,,(0), 

where Sij is a complex number. For the initial state of the system in the computational basis {|00), |01), |10), 
we choose 

m = po|00)(00| +pi|01)(01| +P2|10)(10| 
+(l-Po-Pi-P2)|ll)(ll|+ci2|01)(10| + cf2|10)(01|. 

III. EXACT SOLUTION IN THE MARKOVIAN CASE 

The analytical solution in the basis of eigenvectors |Ai) is given by: 



1 

Pii{t) = 



X1Y2 



/an ai2 ai3 ai4 \ / Pii(O) 

a21 a22 a23 0,24 I I P22{0) 

asi a32 033 034 I I P33(0) 

V a41 042 CI43 044 / \ P44(0) 



where the coefficients are given by: 

an = (X+ + Xre-*^0(^2+ + ^2"e"*''=), 

ai2 = (l-e-*^^)(l-e-*^=)X+y2+, 
ai3 = (1 - e-*^0X+(y2+ + 1^2-6-*''^), 
ai4 = {X+ + e-*^0(l - e-*^=)F+, 

a2i = (l-e-*^^)(l-e~*^^)Xry2-, 



a22 = (Xf + X+e-*^i)(F2- + Y+e-'''-), 
a23 = (Xf + X+e-*^0(l - e-*^=^)F2", 
a24 = (1 - e-'^0^r(^2" + ^2+e-'^^), 
aai = (1 - e-*^0^r(^2+ + >^2"e-*''=), 
a32 = (Xf + X+e-«^)(l - e-*^=)y2+, 

a33 = (Xf + X+e-*^i)(y2+ + F2"e"*''^), 
a34 = (l-e-'^0(l-e"*''=')^r>^2+, 

a42 = (1 - e-«^)X+(y2- + Y+e-'''% 
a43 = (1 - e-*^^)(l - e-*^=)X+F2-, 

044 = (^1+ + Xf e-*^i)(y2" + i"2'^e-*^=). 
Taking into account the initial conditions, the non-vanishing non-diagonal elements are: 

P34(i) = exp f-*2a3 - &±^^ p34(0), 



PAs{t) = P34 = exp ( i2tXz - ^("'^i + -^2) ] p43(o). 



In the solution we have introduced some constants: 



Xi — Xf' + X^ , 



Yi = Y++Y-, 



or 



X=F = 2cos2(^/2)j(i)(±Wi) + 2sin2(0/2)j(2)(±a;i) 
= 2sin2(^/2) J(i)(±a;i) + 2cos^{e /2)J^^\±uJi) 



Xf = (J(i)(±a;,) + J^^\±w,)) + -^=^_(j(i)(±a;,) - /^^(ia;,)) 
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One can easily see that the only steady-state solution possible in this system corresponds to the time moment t = oo: 

( XtY,^ \ 



lim pu{t) 

t—*oo 



X1Y2 



V X+Y,- J 



In the regular basis p^o is: 



lim pzi{t) = 0. 

t^OO 



X1Y2 



/ X{Y^- \ 

c^Xj^Y^ +s^X^Y2' s{Xj;Y^ - X^Y^^) 

s{X{Y^ - X^Y^-) s^Xf K+ +c2X+y2" 

V X+Y+ J 

where c — cos (0/2) and s — sin {0/2). 

In order to quantify the entanglement between the spins we consider the concurrence 
(i ^ 00) it is given by 



Coo = ^Max fo, ^-^\X+Y- - X-Y+\ - ^fx^X+Y^^ 



In the steady-state 



IV. EXACT SOLUTION IN THE POST-MARKOVIAN CASE 

It is a well known fact that positivity is guaranteed only in the case of the Markovian dynamics and in general even 
in the Born-approximation one can find that evolution in no longer positive. Recently Shabani and Lidar [ll|, 
suggested and studied an equation which describes positive and non-Markovian dynamics of the reduced system, so 
called Post-Markovian dynamics 



^ = £ /* dt'k{t') exp {t'C)p{t - t'), 
at Jo 



or 



^ = Ck{t)exp (tC) *p{t). 

Note, that the above dynamics contains a phenomenological memory kernel k{t). 

Solution of the post-Markovian equation can be constructed with the help of the Laplace transform 

C 

sp{s) - p{0) = [k{s) * 

s — L 



Then, the eigenvector-problem for the Lindbladian 



can be solved and one gets the solution 



Cp ~ Xp 



pit) = ^Tr[L,p(t)]i?, =J2^^^it)R^: 
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where 



Hi{t) = Lap 



1 



]m.(o) = C^Wa*.(o). 



s - Xik{s - Aj) 

Particularly, in the case considered in this article the post-Markovian equation takes the following form: 

dp 



dt 



= -i[Hs,p\ + [Ci+ C2) / dt'k{t')cx:p{t'{Ci+ C2))p{t-t'). 



In order to solve the eigenvector problem we find the Jordan decomposition for the Lindbladian 
where 

/ y^/y^- Y+IY2- -I -i\ 

x^/xt -1 X^/X+ -1 

X^Y+IX+Y^- -Y+/Y^ ~X^IX+ 1 
1 1 1 



1 / 



and 



J = diag(0, ~Xi, -Fa, -Xi -Y2). 
In this paper we consider the exponential memory kernel [isj 

which implies that 



The analytical solution for the diagonal elements is 



70 + Ai 



puit) ^SX diag(l,C(J22,t),?(J33,0,C(^44,t)) X S' 



P22(0) 
P33(0) 
V P44(0) 



Taking into account that for the non-diagonal elements the Lindbladian has Jordan form the dynamics of the corre- 
sponding elements are given by ^{Xnon-diag,t) with corresponding value. 



V. RESULTS AND DISCUSSION 

The dynamics of entanglement is analyzed in Figures 1-3. In Figures 1 and 2 the dynamics of the concurrence 
between the two qubits is shown for different coefficients 70 in the memory kernel k{t) (for the Markovian case 
7o = 00). One can see that with decreasing 70 memory effects play a more essential role in the system dynamics and 
practically suppress the oscillations in the concurrence dynamics due to the Hamiltonian dynamics of the system. In 
Figure 3 one can see that increasing the temperature of the baths destroys quantum correlations in the system. From 
Figure 3 and curve (c) on Figure 2 one can see the phenomenon of "sudden death" of entanglement Refs. [IB, 113] ■ 
The steady-state concurrence is analyzed in Figures 4 and 5. The detailed analysis of steady state concurrence for 
this model is given in Ref. 10]. In Figures 4 and 5 we plot the steady-state concurrence for the symmetric (Ae = 0) 
and non-symmetric (Ae 7^ 0) cases as a function the mean temperature (Tm = (?i + T2)/2) and the temperature 
difference (AT — Ti — T2) of the baths. One can see that in the symmetric case (Fig. 4) the maximal value of the 
entanglement corresponds to the thermodynamically equilibrium case (Ti = T2) and in the non-symmetric case (Fig. 
5) the maximum of the quantum correlations reaches in the thermodynamically non-equilibrium case. 

In conclusion, we have found an analytical solution for a simple spin system coupled to bosonic baths at different 
temperatures in Markovian and Post-Markovian cases. We studied the influence of memory effect on the dynamics of 
entanglement. 
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FIG. 1: Dynamics of the concurrence C{t) for the initial reduced density matrix po = |1, 0)(1, 0|. The parameters of the model 

are chosen to be 71 = 72 = 0.001, ei — 2, 62 = 1.1, K = 1, Ti = 0.2, T2 — 0.5 for different coefficients 70 in the memory kernel 
k{t): curve (a) corresponds to the Markovian case k{t) = 5{t); curves (b)-(d) post-Markovian cases; curve (b) 70 = 1; curve (c) 
70 = 0.1; curve (d) 70 = 0.01. 




t Wi ' iJ S i^; 100 12Q 

t 



FIG. 2: Dynamics of the concurrence C{t) for the initial reduced density matrix po = |1, 0){1, 0|. The parameters of the model 
are chosen to be 71 = 72 = 0.001, ei = 2, e2 = 1.1, K = 1, Ti = 1.2, T2 = 1.5 for different coefficients 70 in the memory kernel 
k{t): curve (a) corresponds to Markovian case k{t) = 5{t); curves (b)-(d) post-Markovian cases; curve (b) 70 = 1; curve (c) 
70 = 0.1; curve (d) 70 = 0.01. 
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FIG. 3: Dynamics of the concurrence C{t) for the initial reduced density matrix po = |1, 0)(1, 0|. The parameters of the model 
are chosen to be 71 = 72 = 0.001, ei = 1.5, £2 = 1.1, K = 1, = 0.1 for different temperatures of the "baths": curve (o) 
corresponds to Ti = 0.2, T2 = 0.5; (b) Ti = 1.2, Ta = 1.5; (c) Ti = 2.2, T2 = 2.5. 




FIG. 4: Steady-state concurrence C{Tm, AT) as a function of the mean bath temperature Tm = (Ti + T2)/2 and temperature 
difference AT = Ti — T2 in the symmetric case ei = £2 = 2 with K = 1. 
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